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Aims. We study the stellar population of local infrared galaxies, which contain star-forming galaxies, composite galaxies, LINERs, 
and Seyfert 2s. We also want to find whether infrared luminosity and spectral class have any effect on their stellar populations. 
Methods. The sample galaxies are selected from the main galaxy sample of SDSS-DR4 and then cross-correlated with the IRAS- 
PSCz catalog. We fit our spectra (stellar absorption lines and continua) using the spectral synthesis code STARLIGHT based on the 
templates of simple stellar population and the spectra of star clusters. 
Results. Among the 4 spectral classes, LINERs present the oldest stellar populations, and the other 3 subsamples all present substantial 
young and intermediate age populations and very few old populations. The importance of young populations decreases from star- 
forming, composite, Seyfert 2 to LINER. As for the different infrared luminosity bins, ULIGs & LIGs (log(L 1R /L e ) >11) present 
younger populations than do starbursts and normal galaxies; however, the dominant contributors to mass are old populations in all 
sample galaxies. The fittings also give consistent results by using the spectra of star clusters with different ages and metallicities as 
templates. The dominated populations in star-forming and composite galaxies are those with metallicity Z = 0.2Z o , while LINERs and 
Seyfert 2s are more metal-rich. The normal galaxies are more metal-rich than the ULIGs & LIGs and starbursts for the star-forming 
galaxies within different infrared luminosity bins. 

Additionally, we compare some synthesis results with other parameters obtained from the MPA/JHU catalog. First we find that the 
stellar and nebular extinctions are correlated, and the ionized gas suffers twice as much extinction as stars. Second we confirm that 
D„(4000) is a much better age indicator than HSa- Following the evolution of galaxies, £)„(4000) monotonously varies. Finally we 
investigate some relationships between mean stellar age, mean stellar metallicity, and nebular metallicity for the subsample of star- 
forming galaxies. In star-forming galaxies, the nebular metallicity Z ne t is correlated with the light-weighted mean stellar age ( log t,) L 
in an intermediate strength, and Z„ e t is weakly correlated with the mass-weighted mean stellar metallicity ( Z,) M . 
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1 . Introduction between 3.6 and 160 /-im with better sensitivity and spatial res- 

T . , , . ,. , . r , ■ ■ olution than for IRAS and ISO, and provides new opportuni- 

Infrared galaxies are objects that emit most of their energy in ^ foj . determini the IR propert ies of galaxies in the general 

£j . the infrared (~ 5 - 500 M m) (Sanders & Mirabel 1996). They CQntext of CQsmic evolution ( Le Floc , h et al 2005> zh et 

OS ; are one of the most important discoveries from extragalactic aJ 2QQ7 . p ich et al 2007 and re f er ences therein). The fu- 

observations at mid- and far-infrared wavelengths. The over- ^ Herfjchel wffl alfjo ^ much of Qm knowled about 

whelming majority of the extragalactic objects detected in the infrared laxieg xhereforei in f rare d-detected galaxies are one 

Infrared Astronomical Sate lite (IRAS survey are late-type spi- of ^ mQst interesti objects in the Universe and are related 

if g J ll ™ S ; W ' C : £l ! ipt !. C a l S0 .g alaxief \are rarely detected tQ ^ major star _ formi ss xhe local infrared alaxies 
(Soifer 1987). With the higher sensitivity and spatial resolution . , . eu ■ . ■ • , . •, ... .. 
\ . _ , Tr ^^ „P , J v , . . have the advantage of being studied in details with properties, 
than IRAS, the ISO satellite moves the studies of luminousin- such afj ^ w lationSi since their high . qua ii ty optica i spec - 
frared galaxies to z >0.5 from local (Elbaz et al 1999 2002; ^ d fee obtained . ^ d ^ ^ gk sury ± 
Flores et al. 1999; Liang et al. 2004; Zheng et al. 2004; Hammer S loan Digital Sky Survey (SDSS). 
et al. 2005 and references therein). The cosmic infrared back- 
ground resolved by ISO mid-infrared camera (ISOCAM) shows Understanding the overall stellar population in galaxies 
that the co-moving density of infrared light due to the luminous 18 a cmcial tGo1 for unveiling the star formation and evolu- 
IR galaxies was more than 40-70 times higher at z ~ 1 than tlon of galaxies. As spectra of galaxies record information 
today (Elbaz et al. 2002). And this evolution is mainly driven about a g e and metallicity distributions of its stellar populations, 
by the luminous infrared starburst galaxies, in which the star man Y works have been generated on stellar population analy- 
formation rate runs up to IQQM^yr 1 (Flores et al. 1999). The Sls through spectra of galaxies, using the approach called stellar 
successful launch of the Spitzer Space Telescope opened a new population synthesis. The first attempts were the empirical pop- 
exciting window on the deep infrared universe. Spitzer operates ulation synthesis (Faber 1972), in which the observed spectrum 

of a galaxy is reproduced by a combination of spectra of individ- 

* email: chenxy@bao.ac.cn ual stars or star clusters of different ages and metallicities from 

** email: ycliang@bao.ac.cn a library. More recent models have been based on the evolution- 
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ary population synthesis technique (Tinsley 1978; Bruzual 1983; 
Chariot et al. 1996; Bruzual & Chariot 2003), in which the main 
adjustable parameters could be the stellar initial mass function, 
the star formation rate, and the chemical enrichment. The repro- 
duced simple stellar populations (SSP) with different ages and 
metallicities from the evolutionary population synthesis model 
could be used to analyze the stellar populations of the galaxies. 

Some research has taken on analysis of stellar populations 
in star-forming galaxies and starburst galaxies from the optical 
band observations, most of which conclude that H II galaxies are 
age-composite stellar systems (Schmitt et al. 1996; Raimann et 
al. 2000; Cid Fernandes et al. 2003; Kong et al. 2003; Westera 
et al. 2004). Hammer et al. (2001) used the spectra of star clus- 
ters to analyze the stellar populations of distant luminous blue 
compact galaxies at redshift 0.4 < z < 1 and find that they have 
included some old populations. 

Moreover, some works have been performed on ultraviolet 
(UV) and near infrared bands. Based on International Ultraviolet 
Explore (IUE) UV spectra, Bonatto et al. (2000) discuss the stel- 
lar population of galaxies with enhanced star formation, which 
are at distances corresponding to radial velocities in the range of 
5,000-16,000 km s _1 . They find that young stellar populations 
(age < 500 Myr) are the main contributors except for galaxies 
with a red spectrum arising from the intermediate (age ~ 1-2 
Gyr) and old-age populations. Riffel et al. (2008) have studied 
the central stellar populations of four starburst galaxies through 
near-infrared spectroscopy and find that the dominate population 
was a lGyr old with solar metallicity component in all galaxies. 

Furthermore, some research was applied to AGNs. Cid 
Fernandes et al. present a series of studies that discussed 
the stellar populations of low-luminosity active galactic nu- 
clei (LLAGN) (Cid Fernandes et al. 2004a, 2005a; Gonzalez 
Delgado et al. 2004). They find that young stars contribute 
very little to the optical continuum, while intermediate-age stars 
contribute significantly; and most of the strong-[OI] LLAGNs 
predominantly have old stellar population. Besides that, some 
works have been dedicated to the stellar populations of Seyfert 
2s, and they confirm that most Seyfert 2s contain a young hot 
stellar population (Schmitt et al. 1999; Storchi-Bergmann et al. 
2000; Gonzalez Delgado et al. 2001; Joguet et al. 2001; Cid 
Fernandes et al. 2004b). Boisson et al. (2000) performed stellar 
population synthesis on a sample of 12 galaxies (Boisson et al. 
2004 studied 5 more) with different levels of activity, and they 
find different populations for different activity types: the star- 
burst galaxies present the youngest populations of the sample; 
the Seyfert 2s also have young stellar population, less intense 
than in the starbursts, but metal rich; while LINERs show the 
oldest populations, metal rich, with little star formation still go- 
ing on. But their sample was small, so a much larger sample from 
the modern digital sky survey is needed to verify these results. 

To our knowledge, no work has been specifically dedi- 
cated to the stellar population analysis of local infrared-selected 
galaxies. In this paper, we utilize the data from cross iden- 
tifying between SDSS and IRAS to study the stellar popula- 
tion of local (with median value of redshift ~ 0.03) infrared- 
selected galaxies by using a stellar population synthesis pro- 
gram called STARLIGHT (Cid Fernandes et al. 2005b). With 
the high-quality SDSS optical spectra and the IRAS observa- 
tions, this could be the largest local sample selected from in- 
frared to date for studying their stellar populations from the full 
spectrum fitting pixel by pixel. Moreover, as our sample contains 
diverse spectral classes (star-forming galaxies, composite galax- 
ies, LINERs, and Seyfert 2s) and spans a broad range of infrared 



luminosity, we can compare the properties among these different 
classes. 

This paper is organized as follows. We present the sample 
selection and the classifications in Sect. 2. In Sect. 3 we describe 
the methods of the spectral synthesis and the results of fitting 
the spectra of galaxies in different classes and different bins of 
infrared luminosity. In Sect. 4, we present the relations of some 
property parameters from emission lines, continua of galax- 
ies, and some from the spectral synthesis fittings on the spec- 
tra. Summary and conclusions are given in Sect. 5. Throughout 
this paper we assume the following cosmological parameters: 
Q. M = 0.3, Q. A = 0.7, and H =70 kms" 1 Mpc" 1 . 



2. Sample selection and classification 

Our sample was selected from the main galaxy sample of the 
Sloan Digital Sky Survey (SDSS) DR4 database and then was 
cross-correlated with the IRAS Point Source Catalog (PSCz). 

The SDSS-DR4 includes over 800,000 fiber spectra of ob- 
jects selected based on the five bands (u,g, r, i, and z) imaging 
(Fukugita et al. 1996; Smith et al. 2002; Strauss et al. 2002). 
Each fiber is 3" in diameter (Bernardi et al. 2003; Greene et 
al. 2005, 2006; Adelman-McCarthy et al. 2006), corresponding 
to ~ 1.8 kpc at z = 0.03 (the median redshift of our sample, 
see below and Fig.Q]). The observed-frame spectral wavelength 
range is 3800-9200 A. The instrumental resolution of the spec- 
tra is A/AA * 1800 equivalent to Gaussian cr; nst ss 70kms _1 
(Heckman et al. 2004; Greene et al. 2005, 2006). 

In particular, we selected all galaxies with redshift less than 
0.5 and redshift confidence greater than 0.9 from SDSS DR4. 
After that these objects were cross-identified with PSCz using a 
5" matching radius. The IRAS PSCz survey is a redshift survey 
of about 15,000 galaxies detected, after supplements and cor- 
rections (Saunders et al. 2000), containing 18,351 objects with 
fy) > 0.6/y, covering 84% of the full sky. The optical positions 
of PSCz objects were used to assure the reliability of the identi- 
fication. As the PSCz objects mostly have redshift less than 0. 1 
(Saunders et al. 2000), the redshift distribution of our sample 
objects is almost truncated after a redshift of 0. 1 . 

After cross-identification, 1097 target objects were obtained. 
To improve the accuracy of the classification of our sample ob- 
jects, we excluded 34 objects without emission-line flux mea- 
sured in the MPA/JHU catalog. Meanwhile we excluded another 
65 objects with low emission-line signal-to-noise ratio (S/N) ac- 
cording to their H/3, Ha, and [NII]/16583 lines detected at lower 
than 5cr, and [OIII]/15007 line detected at lower than 3cr. After 
this cut, 998 sources are left. 

To ensure the spectra were taken at the centers of the galaxies 
(within a radius of 2" from the center), we examined the posi- 
tions of fibers for all the sample objects. We found that there 
were 49 galaxies with the fibers not located at the centers. After 
removing them, we have 949 objects left in total. 

To do the spectral synthesis fittings, 100 objects were further 
ejected because their mask spectra were problematic, in which 
almost half of the spectra in wavelength were masked (see more 
details in Sect. |3.1| ). Finally, 849 objects were left as our working 
sample for a stellar population analysis. 

The redshifts of the whole sample of 849 galaxies range from 
0.002 to 0.24, with a median of 0.03, which is presented in Fig.Q] 
(left column). In this work, we used two methods (by emission- 
line ratios and by infrared luminosity) to classify the sample 
galaxies into several subgroups as shown below. 
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2.1. Classification by emission-line ratios 

The first method to classify the samples is using an emission-line 
diagnostic diagram. Namely, we classified the 849 sources as 
star-forming galaxies, composite galaxies, LINERs, and Seyfert 
2s according to the emission-line flux ratio diagnostic diagram 
(Fig© Baldwin et al. 1981, BPT; Veilleux & Osterbrock 1987; 
Kauffmann et al. 2003a; Kewley et al. 2001). First, we sepa- 
rated star-forming galaxies, composite galaxies, and AGNs by 
Eq. © (Kewley et al. 2001, Fig. Ell and Eq. © (Kauffmann 
et al. 2003a, Fig. ©. The composite galaxies refer to objects 
whose spectra contain significant contributions from both AGN 
and star formation (Brinchmann et al. 2004, Kewley et al. 2006). 
Then we divided AGNs into LINERs and Seyfert 2s by Eq. © 
(Shuder et al. 1981). The fluxes of these emission lines were de- 
rived from the MPA/JHU catalog. We finally obtained 419 star- 
forming galaxies, 326 composite galaxies, 69 LINERs, and 35 
Seyfert 2s (Table©: 
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2.2. Classification by infrared luminosity 



Fig. 1. The distributions of redshift (left) and infrared luminosity 
1 .3, (2) (right) f° r our sample of 849 infrared-selected galaxies. 
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We calculated the far-infrared luminosity Lfir of the sample 
galaxies from the nominal flux densities, / v (60 //m) and / v (100 
/mi) in the IRAS catalogs (Helou et al. 1988; Sanders & Mirabel 
1996), and the redshifts given by SDSS, and then we converted 
it to the total infrared luminosity Lir (1-1000 /jm, Calzetti et al. 
2000; Cao et al. 2006; Wang et al. 2006; Wang 2008; Chen et al. 
2008) according to the following formulas: 



F FlR = 1.26 x 10- 14 {2.58/ 60 + /iooHWnT 2 ], 
Lfir = 4nDlF F i R [L Q ], 
L m (l ~ 1000/i m ) = 1.75LHR. 
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The infrared luminosity of the whole sample of 849 galaxies 
spans a range of log(LiR/L ) = 8.65 - 12.53, with a median 
of 10.77, and we show the distribution of infrared luminosity in 
Fig- [U (right column). 

The second classification method was in accordance with 
their infrared luminosities. Namely, we divided our 849 ob- 
jects into 3 subsamples according to their infrared luminosity 
(Elbaz et al. 2002). As a result, there are 299 ULIGs and LIGs 
(L IR /L e > 10 11 ), 451 starbursts (10 10 < L IR /L < 10 11 ), and 99 
normal galaxies (L/«/L < 10 10 ) in our sample. But the num- 
ber of the ULIGs is very small (there are only 20 ULIGs in 849 
objects). 

We divided each of the spectral classes from emission-line 
ratios (i.e. star-forming galaxies, composite galaxies, LINERs, 
and Seyfert 2s) further into 3 infrared luminosity bins as above. 
The related numbers of the objects are given in Table [1] 

3. Spectral synthesis 

In this section, we fit the spectral absorptions and continua of 
the sample galaxies to study their stellar populations by using 



-0.5 

log([NII]\6583/Ha) 



Fig. 2. The emission-line flux ratio diagnostic diagram (i.e. BPT, 
Baldwin et al. 1981) for our full sample (849), in which we plot 
the ratio of [OIII]/H(3 versus the ratio of [Nil]/ Ha . The dotted 
curve defined by Kauffmann et al. (2003a) and the dashed curve 
defined by Kewley et al. (2001) show the separation among star- 
forming galaxies, composite galaxies, and AGNs. The solid line 
defined by Shuder et al. (1981) shows the separation between 
LINERs and Seyfert 2s. 

the software STARLIGHTfJ (Cid Fernandes et al. 2005b, see de- 
tails below). We downloaded the spectra in a one-dimensional 
FITS formaQ, and they have been sky-subtracted, then the tel- 
luric absorption bands were removed before the wavelength and 
spectrophotometric were calibrated (Stoughton et al. 2002). We 
corrected the foreground Galactic extinction using the redden- 
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Table 1. Number of each subsample. 
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class 


star-forming 
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LINER 


Seyfert 2 




total number 


419 


326 


69 


35 


849 


ULIG & LIG 


100 


158 


26 


15 


299 


starburst 


250 


152 


32 


17 


451 


normal 


69 


16 


11 


3 


99 



ing maps of Schlegel et al. (1998) and then shifted the spectra 
to the rest frame. To obtain the general view of the stellar pop- 
ulations of galaxies in different spectral classes and to improve 
the S/N of the spectra, we combined the spectra of all the objects 
to be one spectrum in each of the spectral classes: star-forming 
galaxies, composite galaxies, LINERs, and Seyfert 2s. Likewise, 
to check the effect of infrared luminosity on stellar populations 
of galaxies, we took star-forming galaxies as representative and 
combined their spectra in each of the three infrared luminosity 
bins. We also did the same stellar population synthesis fittings on 
the individual spectrum of the 849 sample galaxies. The results 
will be discussed in Sect.4. 

3.1. The method 

The software STARLIGHT was developed by Cid Fernandes 
and his colleagues to synthesize the stellar absorptions and con- 
tinua of the sample galaxies (Cid Fernandes et al. 2005b, 2007; 
Mateus et al. 2006, Asari et al. 2007). It is a program to fit an 
observed spectrum 0,\ with a model M,\ that adds up to N„ SSPs 
with different ages and metallicities from the evolutionary syn- 
thesis models of Bruzual & Chariot (2003) (hereafter BC03). 
The line-of-sight stellar motions are modeled by a Gaussian dis- 
tribution centered at velocity v» and broadened by cr„. The fit 
is carried out with the Metropolis scheme (Cid Fernandes et al. 
2001), which searches for the minimum^- 2 = T,x[(Ox -Mx)u>x] 2 , 
where is the error in Ox except for masked points. Pixels 
that are more than 3<x away from the rms Ox - Mx are given zero 
weight by the parameter 'clip' . 

To extract the SSP from BC03, we used the "Padova 1994" 
tracks (Alongi et al. 1993; Bressan et al. 1993; Fagotto et al. 
1994a, b; Girardi et al. 1996) and a Chabrier (2003) initial mass 
function (IMF). Extinction was modeled as rising from the fore- 
ground dust, with the reddening law of Calzetti et al. (1994, 
CAL hereafter) with R v = 3.1. As mentioned by the manual 
of STARLIGHT that the reddening law of Cardelli et al. (1989, 
CCM hereafter) was not a good choice for ULIGs, we used CAL 
law. Moreover, we used a base of 45 SSPs, which included 15 
different ages from 1 Myr to 13 Gyr (i.e. 1, 3, 5, 10, 25, 40, 
100, 280, 640, 900 Myr & 1.4, 2.5, 5, 11, 13 Gyr) at each of the 
three metallicities: 0.2, 1, and 2.5 Z . All bases were normalized 
at /lo = 4020A. Additionally, a power law stands for nonstellar 
component (Koski et al. 1978) F Y a V , was added when we fit 
the spectra of AGNs (both LINERs and Seyfert 2s). We set the 
power-law index as a = —1.5 (Cid Fernandes et al. 2004b). We 
will first describe the method for spectra synthesis. All spectra 
were sampled again in steps of A A — 1 A from 3700 to 8000 A, 
and normalized by the median flux in the 4010 to 4060 A region. 
The S/N in the normalization window spans a range of 1 . 1 - 86 . 1 , 
with median value of 15.6. Then we compiled the mask-files 
for our sample. Besides emission line and sky line, we also ex- 
cluded points that have nonzero mask spectraj in SDSS, which 
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are bad pixels or other artifacts. We furthermore excluded four 
other windows: 5870 - 5905 A, to leap the JVa £Ui5890, 5896 
doublet, which is from the interstellar medium; 6845 - 6945A 
and 7550 - 7725 A, which are listed by BC03 as their bugs due to 
problems in these ranges in the STELIB library (Le Borgne et al. 
2003); 7165 - 7210A, which shows a systematic broad residual 
in emission as mentioned by Mateus et al. (2006). In addition, 
there were 100 objects (56 star-forming galaxies, 34 composite 
galaxies, 6 LINERs, and 4 Seyfert 2s) whose mask spectra were 
problematic, in which almost half of the spectra were masked. 
We do not fit these 100 spectra as mentioned in Sect. [2] 

In the outputs of STARLIGHT, one of the most important pa- 
rameters to present stellar population is the population vector x. 
The component xj(J = 1, ...,N*) represents the fractional contri- 
bution of the SSP with age tj and metallicity Zj to the model flux 
at the normalization wavelength Aq = 4020A. Equivalently, an- 
other important parameter, the mass fraction /ij, has the similar 
meaning. We analyze the obtained values of these parameters. 

3.2. The results of fitting the combined spectra according to 
different spectral classes 

As discussed in Sect. 12.11 our sample can be divided into four 
subsamples by a diagnostic diagram of their emission-line ra- 
tios, i.e., star-forming galaxies, composite galaxies, LINERs, 
and Seyfert 2s. To clearly show their properties and to examine 
whether there is any difference in the stellar populations of dif- 
ferent spectral classes, we combine all the spectra of the sample 
galaxies in each subsample, then we synthesize these combined 
spectra with improved S/N. 

Figure [3] shows the spectral fitting results for our combined 
spectra of 419 star-forming galaxies (the four panels at top-left), 
326 composite galaxies (the four panels at top-right), 69 LINERs 
(the four panels at bottom-left), and 35 Seyfert 2s (the four pan- 
els at bottom-right). 

For each of the subgroups, we plot the observed spectra, the 
synthesis spectra, the error spectra, and the residual spectra (the 
pure emission-line spectra) in the left two columns, and the light 
and mass fractions as functions of ages of each SSPs in the right 
columns. We can see that the error spectra are all very weak 
except the area around Hp, and the emission lines, especially 
the weak ones, appear much more clearly in the residual spectra. 

Additionally, at the top right of each subgroup, we also list 
some other parameters, such as;^ 2 , i.e. the reduced^ 2 , the mean 
relative difference between synthesis and observed spectra A^; 
the S/N in the region of 4730-4780 A; the V-band stellar extinc- 
tion; the velocity and the velocity dispersion. Again we confirm 
that our fit is good enough except for the area around Hp. As 
to the bad fit around Hp, which means that there is always a 
trough in the continuum around Hp, Asari et al. (2007) suggests 
that this could relate to calibrations in the STELIB library in this 
spectral range. The good fits around Hp by using the spectra of 
star clusters instead presented in Sect. |3.6| confirm this. However, 
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Fig. 3. The spectral synthesis results of the combined spectra of each spectral class by using 45 SSPs from BC03. The subsamples 
are star-forming galaxies (top-left four panels), composite galaxies (top-right four panels), LINERs (bottom-left four panels), and 
Seyfert 2s (bottom-right four panels). Top left: comparison of synthesis spectrum (red line) with the observed spectrum (black line), 
and green line shows the error spectrum; bottom left: the residual spectrum, green lines represent mask regions given by SDSS flag, 
magenta lines represent mask regions given by ourselves, and red crosses represent clip points judged by STARLIGHT code; top 
right: fraction of light at A® associated to each of the 15 ages of SSPs used in our fits. The red points in this panel of Seyfert 2s and 
LINERs represent the fraction of power-law; bottom right: fraction of mass (log-coordinate) as a function of the 15 ages of SSPs 
used in our fits. In the top right corner of each group, we list some parameters explained in text (please see the online color version 
for more details). 



this problem around H/3 will not affect our general results much 
for the stellar populations of the galaxies. 

To obtain a general view of the stellar populations of the 
sample galaxies, we arrange their stellar populations into three 
age bins: old populations with age > 1 x I0 l0 yr, intermediate 
age populations with age between 6.4 x 10 s yr and 5 x 10 9 yr, 
and young populations with age < 5 x 10 s yr. Then we list the 
contributions of each age bin to the 4 combined spectra in per- 
centage in Table [2] The corresponding fractions for different in- 
frared luminosity bins are also given, and we show the fractions 
of 3 metallicity bins. We should notice that these numbers are the 
percent fractions, which is a bit different from the Xj component 
output by STARLIGHT directly. Combining Tableland the top 
right panels in each subgroup of Fig. [3] we find that star-forming 
galaxies, composite galaxies, and Seyfert 2s contain substantial 
young and intermediate age populations, while old populations 
give very few contributions. As to LINERs, the dominant contri- 
butions come from old and intermediate age populations, while 
there are almost no young components. Therefore we conclude 
that LINERs present the oldest stellar populations among four 
subsamples, and the star-forming galaxies have more younger 



populations than others. Besides that, we find the significance 
of young (< 5 x 10 8 ) populations decreases from star-forming, 
composite, Seyfert 2, to LINER. This gives credibility to our re- 
sults, as Seyfert 2s and LINERs are well known to have their 
activities dominated by AGNs (or possibly shocks for LINERs), 
while starbursts are, in principle, with young stars. These results 
are consistent with Boisson et al. (2000, 2004), Kauffman et al. 
(2003a), Gonzalez Delgado et al. (2004), Cid Fernandes et al. 
(2008), and Stasiriaka et al. (2008). 

Another phenomenon we noticed is that the dominant con- 
tributors to mass are all old stellar populations in four subsam- 
ples, and we speculate that it is due to the much quicker evo- 
lution of young massive stars than old ones. As mentioned in 
Sect. 13.11 we added a featureless continuum (power-law) in fit- 
ting spectra of LINERs and Seyfert 2s, and this component con- 
tributes almost 10% to the light of galaxies (Tabled. For another 
test, we removed the power-law component of LINERs, making 
their young populations increase about 8%. On the other hand, 
if we add a power law to composite galaxies, then the power law 
contributes little (about 3%). There are some changes, but the 
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dominant stellar populations are unchanged, which confirms our 
results. 

The templates that we used in doing spectral synthesis are 
45 in total, including 15 SSPs at each of the three metallicity 
grids Z = O.2Z , Z , 2.5Z . This could let us check the metallic- 
ity effect further, although the degeneracy of metallicity and age 
is popularly known. Then we obtained the contributions of each 
of the three metallicity grids by summing their contributions at 
different ages, which are given in Table[2] It shows that the star- 
forming and composite galaxies are dominated by the popula- 
tions of Z = O.2Z , meaning that most of their populations have 
metallicities even lower than Z , and LINERs and Seyfert 2s are 
more metal-rich. The populations of Z — Z significantly con- 
tribute to LINERs, which is consistent with the suggestion of 
Boisson et al. (2000, 2004), who comment that LINERs show 
the oldest populations in the sample, metal-rich, with little star 
formation still going on. The Seyfert 2s show a significant con- 
tribution from the SSPs with Z = 0.2, 2.5Z , which could be 
consistent with Cid Fernandes et al. (2004b), who suggest that 
the star formation history of Seyfert 2 nuclei is remarkably het- 
erogeneous: young starbursts, intermediate-age, and old stellar 
populations all appear in significant and widely varying propor- 
tions. 

3.3. The results of fitting the combined spectra according to 
infrared luminosity 

To investigate whether the stellar populations of galaxies corre- 
late at all with the infrared luminosity and to avoid the effects 
by spectral class, we specifically combine the spectra of star- 
forming galaxies in three different infrared luminosity bins (as 
mentioned in Sect. 12.2b . and fit their spectra. We list the spectral 
fitting results in Fig. [4] 

Similar to Fig. [3] we get very weak error spectra and 
good residual spectra in left columns except areas around H/3. 
Moreover, in Table [2] we also list the contributions of each age 
bin to the combined spectra. Combining Table [2] and the light 
fraction distributions in Fig. [4] we can see that all subsamples 
contain both remarkable young and intermediate age popula- 
tions and a few old populations, while ULIGs & LIGs present 
the most significant fraction of young populations among three 
subsamples. The population with intermediate age is the domi- 
nant one in the starbursts and normals. As for the stellar popula- 
tions in each grid of metallicities, the normal galaxies are more 
metal-rich than ULIGs & LIGs and starbursts. The last two are 
dominated by the Z = O.2Z components. 

3.4. The spectral synthesis for combined spectra by using 
SSPs of Z B and0.2Z e 

To avoid the age-metallicity degeneracy and to further check the 
robustness of our estimations about the stellar populations of the 
sample galaxies in different age bins, we only adopt the 15 SSPs 
at each of Z = Z and Z = O.2Z to re-do all the fittings on 
the 7 combined spectra studied above. The results are given in 
Tables [3] and |4j respectively. The way we arrange age bins are 
the same as before. The corresponding fractions for different 
infrared luminosity bins are also given. The dominant contri- 
butions have no obvious changes except the 'old' populations 
increase (and the intermediate ones decrease) by some fractions, 
~5-15% in the Z = Z case, and ~ 10-30% in the Z = O.2Z case. 
The corresponding changing fractions are about 1-7% and 2- 
20% for star-forming galaxies in infrared luminosity bins. These 



could be understood simply because the old populations should 
increase their contributions to supply the lack of metal-rich com- 
ponents. Therefore, we believe that our estimates about the stel- 
lar populations of the sample galaxies are robust. 

3.5. The spectral synthesis for typical spectra by using SSPs 

ofZ e 

Up to now, we have been working on the combined spectra of 
galaxies in the related classes. But it would be important to indi- 
cate the variation in our results when a single spectrum instead 
of a synthetic average spectral energy distribution (SED) is used. 
Thus, we took one typical SED in each of the 7 classes and re-do 
all the fits. This time we adopted 15 SSPs with Z = Z following 
Sect. 3.4. We did not expect much differences from the results 
obtained above. The results are given in Table [5] and Figs. |5] [6] 
It is true there are no obvious differences from the results given 
in Table [3] i.e., the case of Z = Z for the combined spectra. 
Also there is not much difference between these and the results 
in Tables [2] and [4] Again, these could confirm our robust analy- 
ses of stellar populations of the sample galaxies. 

3.6. The spectral synthesis by using spectra of star clusters 

In this section, we use 15 spectra of star clusters with different 
ages and metallicities given in Bica & Alloin (1986a,b) to re-fit 
the combined spectra of the subgroup galaxies. The resolutions 
of the two sets of spectra have been matched. The results are 
given in Figs. [7] and [8] for the different spectral classes and in- 
frared luminosity bins. We can see that after changing our bases 
into spectra of the star cluster, the area around H/3 obtains a much 
better fit. 

We then arranged all bases into 4 age-bins: old populations 
> 10 10 yr, young populations < 5 x 10 8 ;yr, intermediate age pop- 
ulations 5 X 10 8 - 10 I( V"i an d H II. We list the stellar population 
fractions of each age bin in all subsamples in Table|6] The same, 
these numbers are the percent fractions. We can see that LINERs 
present an overwhelming fraction from old populations, the star- 
forming galaxies have younger populations than the others, and 
the significance of young (< 5 x 10 8 ) populations decreases from 
star-forming, composite, Seyfert 2, to LINER. These are consis- 
tent with the results obtained by using the BC03's SSP As to the 
star-forming galaxies in different L/« bins, the results are consis- 
tent with those from the BC03's SSP as well. 

In addition, we also used the SSPs of Vazdekis/Miles 
(Vazdekis 1999, Sanchez-Blazquez et al. 2006, Koleva et al. 
2008) instead of those of BC03 to re-fit the spectra, and the re- 
sults are consistent with those from the BC03's SSP. 

4. The relations of some property parameters from 
emission lines, absorption lines, and continua 

In this section, we investigate the relations of some property pa- 
rameters obtained from spectral synthesis and the emission lines 
or continua measurements. 

The outputs of STARLIGHT are divers, and we discuss some 
of them in this paper, i.e., the V-band stellar extinction A* the 
light-weighted mean stellar age ( log f«) i (Eq. [7J, and the mass- 
weighted mean stellar metallicity (Z„) m (Eq.0: 

N 

<logt*>L = 2xjlogt*j, (7) 
j=l 
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Fig. 4. Spectral fittings for the combined spectra of star-forming galaxies in each infrared luminosity bin by using 45 SSPs from 
BC03 . The left four panels are for ULIGs & LIGs ( log{Lm / L ) > 1 1 ), the middle four panels are for starbursts ( 1 < log{LiR / L©) < 
1 1), and the right four panels are for normal galaxies (log(LmlL Q ) < 10). All symbols are the same as in Fig. [3] 

Table 2. Stellar populations in each subsample, i. e., the results of fitting the combined spectra by using 45 SSPs from BC03. 
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Table 3. Stellar populations in each subsample, i. e., results of fitting the combined spectra by using 15 SSPs (Z = Z ) from BC03. 
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Table 4. Stellar populations in each subsample, i. e., results of fitting the combined spectra by using 15 SSPs (Z = O.2Z ) from 
BC03. 
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Notes: All columns have the same content as Table[3] 

Table 5. Stellar populations in each of the 7 typical spectra, i.e., results of fitting typical spectra by using 15 SSPs (Z = Z Q ) from 
BC03. 
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Notes: All columns have the same content as Table[3] 
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Fig. 5. Spectral synthesis of 4 typical spectra of 4 spectral classes by using 15 SSPs (Z = Z ) from BC03 as bases: top-left four 
panels are for star-forming galaxy (PID-MJD-FID: 0310-51990-579), top-right four panels are for composite galaxy (0499-51988- 
631), bottom-left four panels are for LINER (0845-52381-633), and bottom-right four panels are for Seyfert 2 (0933-52642-185). 
All symbols are the same as in Fig. [3] 
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Fig. 6. Spectral synthesis of 3 typical spectra in each Ljr bin of star-forming galaxies by using 15 SSPs (Z = Z ) from BC03 as bases: 
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and the right four panels are for the normal galaxy (0586-52023-131). All symbols are same as Fig. [3] 



N, 

<Z t > M =J^J Z *J' (8) 
j = l 

where the subscript L or M denotes a light-weighted or a mass- 
weighted average. 



The MPA/JHU collaboration has put the measurements of 
emission lines and some property parameters for the SDSS 
galaxies on the MPA SDSS websitsfl (Kauffmann et al. 2003b; 
Brinchmann et al. 2004; Tremonti et al. 2004 etc.). These values 



http://www.mpa-garching.mpg.de/SDSS/ 
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were obtained from the stellar-feature subtracted spectra with the 4. 1 . Relations between stellar extinction and nebular 
spectral population synthesis code of Bruzual & Chariot (2003). extinction 
We used these emission-line flux measurements to obtain some 
of their property parameters, or we adopted some parameters 
from their catalog, e.g., D„(4000) and H8a- 

We estimated the nebular extinction inside the galaxies from the 
Balmer emission-lines ratios Ha/Hfi, assuming case B recombi- 
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Table 6. Stellar populations in each subsample, i. e., results of fitting the combined spectra by using spectra of 15 star clusters as 
bases. 
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Fig. 9. Comparison between stellar and nebular extinctions for star-forming galaxies (left panel), composite galaxies (middle panel), 
LINERs, and Seyfert 2s (right two panels). The red triangles represent ULIGs & LIGs (i.e. L IR /L > 10 11 ), blue open circles 
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sample galaxies there. 



nation, with a density of 100 cm and a temperature of 10 K 
(Osterbrock 1989) following the relation of 

(^Us = (^0)^10-^)-/^)). (9) 

Using the average interstellar extinction law given by 
Osterbrock (1989), we obtained f(Ha) - f(H/3) = -0.37. We 
then calculated A" v eb according to Eqs. [TOlfTTI (Osterbrock 1989, 
Seaton 1979): 

log(Ha/H^) ins = log(Ha/H/?) obs + c x (-0.37), (10) 

Af = cx3. 1/1.47, (11) 

where c is the extinction coefficient, and we set the intrinsic ra- 
tio (Ha/Hf3)j„ s = 2.86 for star-forming galaxies and composite 
galaxies, and (Ha/H/3) ins = 3.1 for AGNs (Veilleux et al. 1995). 

Comparisons between stellar extinction A* v from our spec- 
tral synthesis and nebular extinction Ay calculated by using 
the emission-line flux measurements given in the MPA/JHU 
database for each subsample are presented in Fig. [9] It shows 
that A* v and A'y b are strongly correlated, and the correlation co- 
efficients are very large. 

The left hand panel of Fig. |9]presents the two Ay values of 
star-forming galaxies. The linear fit including all star-forming 
galaxies yields A* = 0.15 + 0.46 x A'y b , which shows that the 
stellar extinction is almost half of the nebular extinction. This 
confirms the results from other authors (Stasihska et al. 2004; 
Cid Fernandes et al. 2005b; Mateus et al. 2006) and accords with 



studies of nearby star-forming galaxies (Calzetti et al. 1994). To 
check whether the correlation between A* v and A"y b is affected 
by infrared luminosity L/«, we assigned different symbols to the 
3 different Ljr bins, and we also plotted the corresponding linear 
least-square fits of them (ULIGs & LIGs: A ' = 0.19+0.44 xA" v eb , 
starbursts: A^ = 0. 12+0.47 xAy b , and the normal galaxies: A* y = 
0.15 + 0.51 x A"y b ). We can see that these lines are close to the 
solid line (slope ±0.02) except for the long-dashed line. Given 
the distributions of stellar and nebular extinctions in the figure, 
we infer that the ionized gas extinction (nebular extinction) may 
be smaller in normal galaxies than in others, so that we get a 
steeper fit result for them. 

Moreover, in the middle panel of Fig. [9] the linear fit includ- 
ing all composite galaxies yields A* v = 0.12 + 0.44 x A n y b , which 
infers that nebular extinction is almost twice as much as stel- 
lar extinction. This conclusion is similar to that of star-forming 
galaxies. Moreover, we mark different symbols for different bins 
of Ljr and plot the linear fit lines for these subsamples. We find 
that all of the three lines are close to the solid line, which is for 
the total sample. 

Finally, in the right two panels of Fig. [9] we perform the 
relations between stellar extinction and nebular extinction for 
LINERs (left panel) and Seyfert 2s (right panel). In the left panel, 
the linear fit containing all LINERs yields A* v - 0.28 + 0.32 x 
A n y k ', which suggests that the ionized gas exceeds twice as much 
extinction as stars. It is different from the results of star-forming 
and composite galaxies. Furthermore, we plot three other linear 
fit lines for the objects in each of the 3 Lm bins in this panel, 
and we find that the behavior of normal galaxies is similar to 
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are shown in each panel for all the sample galaxies there. All the symbols have the same meanings as in Fig. [9] 



normal galaxies in subsamples of star-forming and composite 
galaxies, while the other two subsamples perform differently. In 
the right hand panel, Seyfert 2s show much weaker relationship 
between stellar and nebular extinctions, and the Spearman-Rank 
Order correlation coefficient is smaller than others. Moreover, 
the linear fit including all 35 Seyfert 2s yields A* v = 0.24 + 0.32 x 
Ay. As the number of data points in Seyfert 2s' sub-sample of 
normal galaxies is extremely poor, we did not try linear fit to 
this subsample. There are very few points in AGNs, especially 
in Seyfert 2s, so we cannot draw definite conclusion. In total, 
the ionized gas tends to suffer excess extinction twice as much 
as stars. 

4.2. D„(4000) and H6 A 

We also acquired two age indicators from MPA/JHU database 
directly: £>„(4000) and H6 A . 

D„(4000): The break at 4000A is the strongest discontinu- 
ity in the optical spectra and is available in our sample. It was 
defined by Bruzual (1983) as the ratio of the average flux den- 
sity in the bands 4050 - 4250 and 3750 - 3950 A; however, we 
used narrower continuum bands here (3850-3950 and 4000-4100 
A), which was introduced by Balogh et al. (1999). In this way 
the index D„(4000) is less sensitive to reddening effects. With 
the increasing ages of the stellar populations of galaxies, their 
D„(4000) values will increase as well, which means the higher 
fraction of older populations. 

H5a'- Strong H6 absorption lines arise in galaxies that expe- 
rienced a burst of star formation that ended ~0. 1-1 Gyr ago. The 
peak occurs once hot O and B stars, which have weak intrin- 
sic absorption, have terminated their evolution. The optical light 
from the galaxies is then dominated by late-B to early-F stars. 
Worthey & Ottaviani (1997) defined an HSa index, which con- 
tains a central bandpass (4083 - 4122A in MPA/JHU database) 
bracketed by two pseudo-continuum band-passes. 

The relationships between age indicators, D„(4000) and 
H5a, and the light-weighted mean age (logt*) l for each sub- 
sample are performed in Fig. \W\ We also show the Spearman- 
rank order correlation coefficient in each panel. In the left 
two panels, we find there is a very strong correlation between 
Z)„(4000) and our light- weighted mean age for star-forming 
galaxies. Another age indicator, H5a, presents a weaker cor- 



relation (lower correlation coefficient) with mean stellar age 
than does D„(4000). The weaker correlation between H6a and 
the mean stellar age may be due to the fact that H5a does 
not monotonously vary following the evolution of the galaxy 
(Kauffmann et al. 2003b). Therefore, we suggest that D„(4000) 
is a good age indicator as the scatter is small and the relationship 
is tight in the left hand panel of Fig.fTOl 

For composite galaxies (middle two panels in Fig. ITOb . 
we draw a similar conclusion to that of star-forming galaxies: 
Z)„(4000) is a good age indicator, as the correlation coefficient 
is high and the diffusion is small. Although H5 A performs better 
than that of star-forming galaxies, there are still some diffusions 
in the figure. 

The right four panels show the results of LINERs and Seyfert 
2s: the left two for LINERs and the right two for Seyfert 2s. 
In the left two we find there are strong correlations between 
D„(4000), Hd A , and ( log t,) L for LINERs, while from the right 
two we find D„(4000) is still a good age indicator for Seyfert 2s. 
In view of so few data points contained in LINERs and Seyfert 
2s, especially Seyfert 2s, we cannot draw any absolute conclu- 
sion. 

Figure [TOl generally shows that the ULIGs and LIGs have 
relatively lower light-weighted mean ages than others, which is 
consistent with our results above (Sect. 3.3). In other words, we 
have tested all individual galaxies, so that we can confirm the 
robustness of the results from fitting the combined spectra. 

4.3. Nebular metallicity and mean stellar age, mean stellar 
metallicity 

In this section, we show the distributions between nebular metal- 
licity Z ne i, and the light-weighted mean stellar age ( log f„) i 
(Eq. |7j, the mass -weighted mean stellar metallicity ( m 
(Eq. |8) provided by STARLIGHT. We only show the results of 
star-forming galaxies for two reasons, one of which is the un- 
certainty of the estimation of Z ne t, for composite galaxies and 
AGNs. Another reason is that we want to compare our results 
with that of Asari et al. (2007), whose sample only includes all 
the star-forming galaxies. 

After correcting all of the useful line fluxes for extinc- 
tion, we quantified the nebular metallicity Z ne b by the oxy- 
gen abundance of the star-forming galaxies. To be consis- 
tent with Asari et al. (2007), the nebular metallicity of star- 
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and the mass-weighted mean stellar metallicity ( Z*) m (right panel) for star-forming galaxies. The Spearman-rank order correlation 
coefficient are shown in each panel for all the sample galaxies there. All the symbols have the same meanings as in Fig. [9] 



forming galaxies was obtained according to Eq. ( TT~2T > where 
O3^ 2 =log([OIII]^5007/[NII] / 16583) (Stasinska et al. 2006; 
Pettini et al. 2004; Liang et al. 2006): 

logZ°*f = log(0/H)/log(0/H) Q = -0. 14-0.25 x0 3 N 2 [Z ],(12) 

where (O/H) =4.9xlO- 4 . 

In Fig. QT| we show the distributions between Z„ e t, and 
(logU)i (the left hand panel), (Z»)m (the right hand panel). 
From the left hand panel, we find there is a correlation with inter- 
mediate strength, and the correlation coefficient is 0.32. It seems 
that the young stellar population makes a relatively dominant 
contribution in low metallicity galaxies. The right hand panel 
shows that there is a weak correlation between the mean stellar 
metallicity and the nebular metallicity, and the correlation co- 
efficient is only 0.17. The chemical enrichment histories of star 
and gas are correlated on the basis of simple chemical evolution 
scenarios, while the weakness and scatter of this relation may 
be due to either the disturbance from interactions/mergers and 
the indefinability of our calculations or the intrinsically different 
evolutionary phases (Asari et al. 2007). These two relations are 
similar to those of the total star-forming galaxies presented in 
the Fig.5 in Asari et al. (2007). 

5. Summary and conclusions 

Based on a large sample (849) of local infrared galaxies se- 
lected from the cross-identification of the SDSS DR4 and IRAS 
PSCz, we fit the spectral absorptions and continua of the sam- 
ple galaxies to study their stellar populations by using the SSPs 
of BC03 and spectra of star clusters as bases and the software 
STARLIGHT. The sample galaxies are divided into four or three 
subgroups by following their emission-line ratios or infrared lu- 
minosity. First, we investigated the effects from spectral class 
(star-forming, composite galaxies, LINERs, and Seyfert 2s) and 
infrared luminosity (ULIGs & LIGs, starbursts, normal galaxies) 
on stellar populations. Second, we fit our sample galaxies with 
SSPs from BC03 and spectra of star clusters. Finally, we com- 
pared the stellar extinction, mean stellar age, and mean stellar 
metallicity derived from spectral synthesis with nebular extinc- 
tion, nebular metallicity, and two age indicators, D„(4000) and 
H5a, obtained from the MPA/JHU catalog. 

As to the contributions to light, we found that LINERs pre- 
sented the oldest stellar populations among 4 spectral classes, 



and the importance of young populations was decreasing from 
star-forming, composite galaxies, Seyfert 2s, to LINERs. The 
dominant populations of star-forming and composite galaxies 
are those with metallicity Z = 0.2Z o , and the LINERs and 
Seyfert 2s are more metal-rich. While ULIGs & LIGs present 
the youngest populations among the star-forming galaxies in 3 
infrared luminosity bins. The normal galaxies are more metal- 
rich than the ULIGs & LIGs and starbursts. However, the domi- 
nant contributors to mass are all old stellar populations, and the 
spectral class and the infrared luminosity have no effect on this 
result. The results of fittings by using spectra of star clusters as 
bases are consistent with those from the SSPs of BC03. 

The stellar and nebular extinctions are correlated in all sam- 
ple galaxies, and the ionized gas suffers twice as much extinction 
as a star. The D„(4000) is a much better age indicator than HSa, 
which does not monotonously vary following the evolution of 
the galaxy. When we compared nebular metallicity with other 
parameters for star-forming galaxies, we found that the mean 
stellar age is correlated with nebular metallicity in an intermedi- 
ate strength, and Z„ e b is weakly correlated with (Z„) m- 
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